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The strong light-matter interaction in graphene over a broad frequency range has opened up a 
plethora of photonics applications of graphene. The goal of this paper is to present the voltage 
tunability of plasmons in gated single- and multi-layer graphene structures. Device concepts for 
plasmonic interconnects and antennas and their performance for THz communication are presented. 

For the first time, the role of gate voltage and the thickness of the gate dielectric on the character¬ 
istics of plasmon propagation in graphene are quantified by accounting for both the interface trap 
capacitance and the quantum capacitance. The gate voltage serves as a powerful knob to tweak the 
carrier concentration and allows building electrically reconhgurable terahertz devices. By optimiz¬ 
ing the gate voltage to maximize the plasmon propagation length in a gated multi-layer graphene 
geometry, we derive simple scaling trends that give intuitive insight into device modeling and design. 


I. Introduction 

The two-dimensional material graphene, which is a 
layer of carbon atoms arranged in a honey-comb lattice, 
exhibits strong light-matter interaction over a broad fre¬ 
quency spectrum from the far infrared to the ultravio- 
lePai. The tunability of the density-of-states in graphene 
along with its excellent transport properties reflected 
in a high carrier mobility provide a path for graphene 
photonic applications such as quantum optics^, photo¬ 
voltaic^, photo-detector^D, and biological sensing 89 . 
The unique optical properties of graphene originate from 
a combination of its two-dimensional nature and gap¬ 
less electronic spectrum. In the optical frequency range, 
the universal dynamical conductivity in graphene 32 is 
(Tmin = 7re 2 /2 h, where e is the elementary charge, and 
h is the Planck’s constant. This provides mono-layer 
graphene with a transparency T = (1 + 2ir(j rn in/c )~ 2 ~ 
1 — 7ra = 0.977, where a = e 2 /(he) is the fine-structure 
constant, h is the reduced Planck’s constant, and c is the 
speed of light. Hence, graphene has been sought after 
as a transparent electrode in solar cell and other flexible 
electronics applications 10 . 

Further, the excitation of plasmons, which are the 
quantized collective oscillations of charged carriers at the 
interface between graphene and dielectric, are responsible 
for several interesting optical attributes of graphene 11 . 
Plasmons lead to the dramatic alteration in optical ab¬ 
sorption features in graphene and can be localized within 
small regions compared to the wavelength of the inci¬ 
dent radiation. In graphene, the propagation length of 
plasmons can be several micrometers and their propa¬ 
gation velocity has a lower bound of Vf/ 2, where vj = 
8 x 10 5 m/s is the Fermi velocity of the Dirac fermions 
in graphene. The long plasmon lifetime and their very 
high propagation velocity make graphene an ideal plat¬ 
form for implementing plasmonic waveguides for on-chip 
communication and ultra-broadband antennas for wire¬ 
less communication. To design photonic devices using 
graphene, analytical descriptions of optical properties in 
graphene prove useful. 


We begin with a discussion of the plasmon dispersion 
relationship for both ungated and gated graphene struc¬ 
tures in Section II. This section details the key differ¬ 
ences in the propagation of surface plasmons in graphene 
upon the application of a gate voltage. In Section III, 
we present the functional relationship between the Fermi 
level (carrier concentration) in graphene and the gate 
voltage using a circuit model to account for trap charges 
due to dangling bonds at the interface between graphene 
and the substrate. In Section IV, we combine the re¬ 
sults from Sections II and III to quantitatively predict 
figures of merit for graphene plasmonic waveguides and 
antennas for wireless THz communication. The paper 
concludes in Section IV with a summary of the results 
presented in the paper and an outlook on the future of 
graphene plasmonics. We also supplement the work with 
three appendices that detail out the specific mathemati¬ 
cal calculations underlying the theoretical foundations of 
the paper. 

II. Plasmon dispersion relation in single- and 
multi-layer graphene 

The plasmon dispersion relationship in graphene is de¬ 
rived from finding the root s of t he real part of the dielec¬ 
tric constant, e Ef , given as 12 13 

e Ef (q, u)=e- v q Ml Ef (q, u) = 0, (1) 

where Ef is the Fermi level in graphene, e is the per¬ 
mittivity of the surrounding media, v q = e 2 /2eo q is the 
Coulomb interaction term, and U Ef (g, cj) is the polariza¬ 
tion function of doped graphene. In the above equation, 
q denotes the wave-vector, and uj is the frequency. Since 
v q > 0, (fTl) has a solution only when ^RU Ef (q,uj) > 0. It 
can be shown that the inequality is satisfied only when 
uj > Vfq. In violation of this condition, one arrives at 
electron-hole continuum in the regime of single-particle 
excitation, which precludes the existence of plasmons 14 . 
Ignoring the non-local effects and in the long wavelength 
limit (q 0) allows approximating the polarization func- 
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tion in doped graphene according to 

We , M = ^A(±) 2 . (2) 

In this case, a simplified plasmon dispersion relation 
in graphene is given ad 16 ^ 


u(o) = 


1 j e 2 \Ef\q ' 
h V 27reoe 


( 3 ) 


Alternatively, the plasmon dispersion relationship can 
be derived by using Maxwell’s equations and applying ap¬ 
propriate boundary conditions for a TM-polarized elec¬ 
tromagnetic wave interacting with graphene embedded in 
a dielectric environment within the non-retarded regime 
where q » uj/c. More details on the derivation of the 
plasmon dispersion relationship for TM-polarized plas¬ 
mon modes in graphene are given in Appendix A. The 
plasmon dispersion relation in graphene is distinctly dif¬ 
ferent from other 2D materials where the dependence of 
v fq is certainly obtained, but the dependence on ^/E^ is 
unique to graphene because of the presence of massless 
Dirac fermions at Brillouin zone edges K and K . Us¬ 
ing the plasmon propagation velocity, v g = duo/dq is 
given as 


Vg(q) = 


1 

h 


' e 2 \E 


71 


87reoeg 


( 4 ) 


Noting that the assumption for the derivation of the dis¬ 
persion relation in Q is that uo > Vfq , it can be shown 
using and @ that the lower bound on the propaga¬ 
tion velocity of graphene plasmons is Vf / 2. 

The propagation length of plasmons in graphene is 
given as 


Lprop — p y (5) 

where the factor T is a phenomenological parameter that 
characterizes the electron scattering rate in graphene 16 . 
An estimate of T can be obtained from the D.C. relax¬ 
ation time, which arises mainly from scattering due to 
intrinsic phonons and the charged impurities of the sub¬ 
strate. In this work, we consider both scattering events, 
which are carrier concentration dependent, to obtain the 
net electron scattering rate in graphene. Details of var¬ 
ious electron scattering times in graphene are noted in 
Appendix B, while the modulation of the carrier con¬ 
centration with an external gate voltage is discussed in 
Section III. 

A. Impact of metal top gate on plasmon dispersion 
in graphene 

When a metal top gate is placed on the graphene sheet 
with a gate dielectric as shown in Fig. [lj the electric field 
distribution is modified such that the electric field has 
only the z— component in the direction perpendicular 


to the graphene sheet as noted in the figure. This is in 
stark contrast to the electric field distribution in ungated 
graphene structure, where electric field has components 
in both z— and x— directions, where the x— direction 
is along the length of the graphene ribbon and coincides 
with the direction of plasmon propagation. Upon appli¬ 
cation of the gate voltage, the plasmon dispersion rela¬ 
tion in graphene is modified according to 


u(q) 



q 

1 + coth qd ’ 


(6) 


where cro is the D.C. conductivity of the graphene sheet 
and is dependent on the Fermi level in graphene, T is the 
electron scattering rate, and d is gate dielectric thick¬ 
ness. The plasmon dispersion relation given above was 
first discussed by Nakayama and co-workers in the con¬ 
text of plasmons in 2D electron gas 17 , where the con¬ 
ductivity in (|6| stands for the conductivity of the 2D 
electron gas. Later, the same dispersion relationship for 
plasmons in gated graphene structures was derived for 
graphene irP^p- 9 and discussed in detail in 20 . 

Using the random-phase approximation in the long 
wavelength limit (q —> 0), the frequency-dependent intra¬ 
band conductivity of mono-layer graphene can be ex¬ 
pressed ad^H 


&intra(^) — 


e 2 cj f°° dE\E\ dfp(E) 

iirh J 0 uj 2 dE 


( 7 ) 


where fo(E) = (1 + exp ((E — E^/ksT)) -1 is the Fermi 
function. Carrying out the integration, cr^ n t ra (^), can be 
written more simply as 


&intra(^) — 


. 2e 2 k B T 
1 1rh 2 (cu + iT) 


In 


2 cosh 


M" 

k B T_ 


( 8 ) 


where we have replaced uj with uj + iT to account for the 
electron scattering rate in graphene. Further, using the 
simplification for the function cosh (x), we can show that 
the intra band conductivity reduces to the form: 


&intra (^) 


. 2e 2 k B T 
irh 2 uj + iT 


Ef 

2k B T 


+ In 


(l + e -E f /k B T 



The D.C. conductivity from the above equation is found 
by setting uj « Y and is given as 




2e 2 k B T 


Ef 

2 k B T 


+ In (l+e~ E f /kBT ^ 


( 10 ) 


In the zero temperature limit, T « 0, which is a good ap¬ 
proximation for heavily-doped graphene (Ef » k B T ), 




e 2 |£/l 

nh 2 T ' 


( 11 ) 


Using this definition of cro in uo(q) can be simplified 
as 


w(<?) = 


i .IfMJ 

h V 7reoe y l + cothgd 


( 12 ) 
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In the limit qd oo, (12) gives the same results as ©• 
In the limit qd « 1, uJ{q) in (12) approaches 


u(q) 


1 

h 



q , M « !)• 


(13) 


As evident from the above equation, in the gated 
graphene structure, plasmon-dispersion relation depends 
linearly on the wave-vector, while in the ungated 
graphene structure the plasmon dispersion relation varies 
as y/q. That is, in gated graphene structures, plasma 
waves exhibit a sound-wave-like dispersion relation. The 
frequency of these plasma waves can be tuned by gate 
voltage within a fairly wide frequency band. 


,v 


d ,l. 

i + + + +-. 


+ + + + 
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FIG. 1: Gated graphene structure with a top gate with a bias 
V g leads to altered plasmon dispersion relation. 

The difference between the dispersion relations in 
gated and ungated graphene structures translates to a 
significant difference in their plasmon propagation veloc¬ 
ity. Using |6]), the plasmon propagation velocity, v gi is 
given as 


Vg 2 1 


a 0 Tq 


eeo (1 + coth (qd)) [qd 1 + cothgd 


1 ^ cschgcr 


• ( 14 ) 


In the limit qd « 1, the plasmon propagation velocity 
approaches the limit 


(JoTd 

eoe 


(15) 


Note that, in this case, the plasmon propagation veloc¬ 
ity is independent of the wave-vector and is a constant 
that depends on the material properties and the geomet¬ 
ric dimension d. This behavior is significantly different 
from that described in where the propagation veloc¬ 
ity scales as ^/g, a characteristic feature of 2D electron 
gas. As shown in Fig. [2j the plasmon propagation ve¬ 
locity for a gated graphene structure converges to that 
for the ungated graphene structure in the limit qd oo. 
In the other extreme when qd « 1, the propagation ve¬ 
locity is simply given as (15), and in this case, the prop¬ 


agation velocity is a constant independent of the wave- 
vector. 


It must also be noted that in the limit qd « 1, 
the propagation velocity of plasmons is lower than the 
propagation velocity of ungated plasmons. For a fixed 
Fermi-level, in this regime, the plasmon propagation ve¬ 
locity scales as \[d, such that the plasmon propagation 
velocity in gated graphene structures with very thin gate 
dielectrics may be significantly lower than the plasmon 
propagation velocity in ungated graphene structres as 
shown in Fig. [3j Despite the lower propagation velocity, 
the advantage of the gated graphene structure is to pro¬ 
vide dynamic tunability of the plasmon characteristics 
in graphene depending on the required circuit function¬ 
ality. This “reconfigurability” option is a fundamental 
advantage of graphene plasmonics over metal-based plas- 
monics. As discussed in Section III, the gate voltage is 
used to tune the carrier concentration and, therefore, the 
Fermi level in the graphene sheet. 



FIG. 2: Plasmon propagation velocity in mono-layer graphene 
normalized to the speed of light versus frequency of operation. 
Three cases are considered: (i) ungated graphene structure 
for which v g is described by Q , (ii) gated grap hene structure 
in the general case for which v g described by (14), and (iii) 
gated graphene structure for which qd « 1 such that v g is 
independent of the wave vector and, therefore, the frequency 
of operation as described by (15). Here, Ef denotes the Fermi 
level in graphene, d is the thickness of the gate dielectric, and 
e is the average dielectric of the surrounding media. 


We would like to note that the plasmon dispersion 
relation derived for mono-layer graphene in the THz 
frequency will be preserved in the case of multi-layer 
graphene as long as the plasma wavelength is much larger 
than the thickness of the sheefP^. However, the carrier 
concentration must include contribution from each layer 
in the multi-layer stack. For perfectly coupled layers 
considered in this work, cj multi = iV Zayer cr mono , where 
(Tmulti (&mono) is the multi- (mono-) layer conductivity 
of graphene, and Ni ayer is the number of layers in the 
multi-layer stack. 
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FIG. 3: Plasmon propagation velocity versus frequency of 
operation for different thicknesses of the gate dielectric, d. 


III. Modulation of Fermi level via gate voltage in 

graphene 


To determine the modulation of carrier concentration, 
n s , in the graphene sheet due to the gate voltage, V g , the 
capacitance voltage network shown in Fig. [4] is used. As 
shown, the equivalent circuit accounts for the capacitance 
due to interface traps created due to dangling bonds at 
the interface between graphene and the substrate. From 
the equivalent circuit diagram, the relationship between 
the Fermi level, Ef , in the graphene sheet and the gate 
voltage, V g , can be expressed as 



e 


C 0 


Cox + C s i + 1/2 Cq 


( V 9 - f'rmn) ■ 


(16) 


Here, C s i is given as the series combination of C c \ and 
C lx , where C c \ is the electrostatic capacitance of the in¬ 
terface traps with the graphene sheet, and C lx is the 
quantum capacitance of the interface traps. C 7l is given 
as e 2 7 i with 71 being the energy-independent density- 
of-states of the interface traps. Typical values of 71 are 
« 5 x 10 12 eV~ 1 cm~ 2 as measured experimentally in 22 
and 23 and also considered by authors in 2 ^ to fit their the¬ 
oretical model to explain the dependence of n s on V g . In 
the above equation, Vmi n is the Dirac point voltage of 
the graphene sheet, and C q = e 2 dn s /dEf is the quantum 
capacitance of graphene. 

To relate the carrier concentration, n s , with the Fermi 
level, Ef , in the graphene sheet, the following Fermi- 
Dirac integrals are used: 


where DOS(E) is the 2D density-of-states in graphene. 
In the presence of electron-hole puddles, the density-of- 
states is more appropriately given ad 25 


DOS(E) 


2 1 - 2(7 dis 

7r(hv f ) 2 y/2ir 


exp 


+E x erf 


E 2 

2cr dis 

E 


Odis 


V2 


(18) 


where is the broadening in eV of the DOS around 
the Dirac point. Using © ~([l 8 |), we show the depen¬ 
dence of E f and n s on the gate voltage in Fig. [5j As 
expected, with an increase in the interface traps DOS, 
the dependence of Ef on the gate voltage is weakened, 
since for the same gate voltage a large amount of capac¬ 
itance comprising of the parallel combination of C s 1 and 
1/2 C q will need to be charged reducing the Fermi level 
in the graphene sheet. An increase in the broadening in 
the DOS around the Dirac point enhances the net sheet 
charge concentration in the graphene sheet. 


IV. Gate-voltage-controlled plasmon propagation 
characteristics in mono- and multi-layer 
graphene ribbons and implications for THz 
wireless communication 


In this section, we quantify the limits of plasmon-based 
on-chip communication by utilizing the gate-voltage tun- 
ability of the plasmon dispersion relationship in both 
mono- and multi-layer graphene structures. We also 
present a preliminary analysis of the frequency response 
of patch antennas with graphene targeted toward several 
interesting communication and radar applications in the 
THz band. Graphene patch antennas have previously 
been discussed inP3. However, prior work focuses on nar¬ 
row graphene ribbons where the Fermi level is only chemi¬ 
cally tunable via doping the graphene sheet. Also, earlier 
works such on graphene plasmonic on-chip intercon¬ 
nect analysis focus only on mono-layer and chemically- 
doped graphene ribbons. Two useful figures of merit are 
discussed in detail in the next sub-sections (i) plasmon 
propagation velocity that determines the latency of high¬ 
speed on-chip interconnects, and (ii) plasmon propaga¬ 
tion length relating the energy per bit for communica¬ 
tion and the radiation efficiency of plasmonic antennas. 
A large propagation length signifies lower ohmic losses in 
graphene and, therefore, superior radiation efficiency of 
the antenna. 


Tls l^eZecl^/) ^hole (-^/) 1 1 


(17a) 


A. Propagation characteristics versus gate voltage 
in 2D graphene 


W'elec 


^hole 


POO 

(Ef) = / 

Jo 

POO 

(Ef) = 

Jo 


DOS(E) 


DOS(E) 


dE 

l + ex p( E k if 1 )’ 

dE 

1 + ex P ’ 


, , In a gated graphene structure, propagation character- 

' ' istics of plasmons are tunable by changing the gate volt¬ 

age. In Fig. [ 6 ] we show plasmon propagation velocity 
versus gate voltage for different values of the number of 
graphene layers and varying dielectric thicknesses. Due 
(17c) to an increase in the 2D conductivity of the graphene 
sheet with an increase in the number of layers in the 
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FIG. 4: Schematic of a graphene sheet on a substrate with dangling bonds. Equivalent circuit model that acounts for interface 
traps due to dangling bonds is also shown. C c i and C 71 are the electrostatic and quantum capacitance, respectively, of the 
interface traps. C q is the quantum capacitance of the graphene sheet. 




FIG. 5: (left) Dependence of Fermi level and (right) sheet carrier concentration on gate voltage. In the left hand side plot, 
different traces correspond to the different values of the density-of-states, 71 of the interface traps. In the right hand side figure, 
different traces correspond to the different values of the broadening in the density-of-states, crais, around the Dirac point in 
graphene. 


multi-layer stack, the propagation velocity is superior for 
Niayer = 5 versus Ni ayer = 1. However, an important 
difference must be noted between the results in Fig. [6] 
compared to the results in Fig. [3] Both figures display 
an opposite dependence on the parameter d, which is the 
gate dielectric thickness. The reason is that in Fig. [6] an 
increase in d leads to a reduction in the Fermi level in 
the graphene sheet because of the reduced gate control. 
This leads to a reduction in the conductivity, cr 0 , which 
in turn reduces the plasmon propagation velocity. In Fig. 

t we had assumed a constant value of Ef in the graphene 
eet and Ef and d were adjusted independently of each 
other. 

We show the dependence of the plasmon propagation 
length on the thickness of the top gate for different gate 
voltages and different number of layers in a gated multi¬ 
layer graphene stack in Fig. [7J As the thickness d in¬ 
creases, the propagation length drops and becomes inde¬ 
pendent of d for very large value of d. This is expected 
since in the limiting case of qd oo, the dependence 
of v g on d vanishes. This makes the plasmon propaga¬ 
tion length independent of d in the limiting case. For the 
same value of relaxation rate, T, propagation length im¬ 
proves with the number of layers in the multi-layer stack. 
It must be noted that while the propagation velocity in- 



Gate voltage (V) 

FIG. 6 : Gate voltage impact on the propagation velocity of 
plasmons in 2D graphene with Ni ayer — 1 (mono-layer) and 
Niayer = 5 (multi-layer) with different values of the gate di¬ 
electric thickness. 


creases monotonically with an increase in the Fermi level 
in graphene, the propagation length of plasmons is a non¬ 
monotonic function of the Fermi level. This is explained 
by observing that an increase in Ef enhances both v g 
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and T such that the propagation length as defined in (J5| 
exhibits a non-monotonic dependence on Ef and, there¬ 
fore, on the gate voltage V g as shown in the inset plot of 
Fig. 0 



FIG. 7: Propagation length versus the gate dielectric thick¬ 
ness. The inset plot shows the dependence of the propagation 
length on the gate voltage. 

As shown in Fig. [8j an optimization of the gate volt¬ 
age to maximize plasmon propagation length reveals two 
noteworthy features. First, the optimal voltage is inde¬ 
pendent of the number of layers in the multi-layer stack 
if the electron relaxation rate is independent of the num¬ 
ber of layers. While the optimal gate voltage to maximize 
the propagation length scales linearly with the value of 
the gate dielectric thickness, the maximum value of the 
propagation length stays constant for a given number of 
layers in the multi-layer stack. This allows us to identify 
a simple scaling law for optimizing the gate voltage to 
obtain maximum plasmon propagation for the graphene 
structure. The scaling law can be stated as 

Vg,opt(d) = Vg,opt(do) + A(d — do), (19) 

where A depends on the relative strengths of the dif¬ 
ferent scattering processes in graphene, and do is the 
reference dielectric thickness. Typical value of A « (70- 
100) mV/nm for a impurity concentration of (5-10) x 10 14 
era -2 . The value of A increases with an increase in the 
impurity concentration. Details of the derivation of the 
scaling law are provided in Appendix C. 

B. Microstrip antennas with graphene plasmon 
resonant cavity 

Microstrip or patch antennas are used in several com¬ 
munication applications owing to their simplicity, ease of 
fabrication, and unidirectional radiation. The advantage 
of graphene to implement high-frequency microstrip an¬ 
tennas is its planar structure and flexibility to be trans¬ 
ferred on several different substrates. This eases the 
integration of graphene within the existing semiconduc¬ 
tor nanotechnologies spanning from biological and chem¬ 



FIG. 8: Gate voltage dependence of the propagation length. 
An optimal gate voltage is identified that maximizes L pr0 p for 
a given number of layers in the multi-layer stack. The value 
of the optimal gate voltage is independent of the number of 
layers. Maximum propagation length is independent of the 
gate voltage. 


ical nanosensor networks to optical interconnects in ad¬ 
vanced multi-core architectures. In this work, we model 
graphene microstrip as a plasmon resonant cavity. As dis¬ 
cussed iiP, a resonant cavity of graphene will necessitate 
the following dimensional constraints to be satisfied 

L » W » h, (20) 


where L and W are the length and the width of the res¬ 
onant cavity, while h is the height of the resonant cavity 
from the substrate. Further, the length of the resonant 
cavity must be selected for a given resonant frequency 
response. For a TM plasmon propagation mode, which 
is the focus of this work, the length L of the graphene 
resonant patch must satisfy 


L 



7r 



( 21 ) 


where m = 1 , 2 ,...., and X spp is the wave-vector of the 
propagating plasmon mode, and q is the plasmon wave- 
vector as discussed in the previous sub-sections. For the 
first TM propagation mode with m = 1, we plot the res¬ 
onant frequency versus the length of the graphene sheet 
for different layers in the graphene multi-layer stack in 
Fig. 0 Both gated and ungated graphene structures 
are considered. As seen, the resonant frequency spans 
from few hundreds of GHz to a few THz depending on 
the length of the patch and the position of the Fermi 
level, which has been adjusted using the gate voltage. 
If one were to utilize chemical doping to achieve the de¬ 
sired Fermi level in the graphene sheet such that the plas¬ 
mon dispersion relation is expressed using for ungated 
graphene structure, the resonant freqency scales as 1 /y/L 
as opposed to scaling as 1 /L for gate-voltage-controlled 
Fermi level. Thus, the resonant frequency for the same 
graphene patch length is lower for gated graphene struc¬ 
ture than for the ungated graphene structure. However, 
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FIG. 9: Resonant frequency as a function of patch length for 
gated and ungated graphene microstrip antennas. 


it must be noted that in either of the two cases, only 
a very small patch length of graphene would suffice for 
THz radiation. 


V. Conclusions 

In this work, we propose device concepts in the tera¬ 
hertz frequency region using plasma oscillations in gated 
graphene structures. In such structures, it is shown that 
sufficiently long plasma waves exhibit a linear (sound¬ 
like) dispersion. Since graphene supports mass less Dirac 
fermions, the plasma waves in graphene can fall within 
the terahertz band even though their wavelength is suf¬ 
ficiently long. By modeling graphene as a resonant cav¬ 
ity or a voltage-controlled waveguide, different voltage- 
tunable terahertz devices can be materialized. We char¬ 
acterize the performance of on-chip interconnects and 
terahertz antennas built using mono- and multi-layer 
graphene patch on a substrate. Accounting for both 
intrinsic acoustic phonons and charged impurities scat¬ 
tering, we derive a simple scaling law between the ap¬ 
plied gate voltage (“reconhgurable” knob) and the thick¬ 
ness of the gate dielectric (“fixed” upon manufacturing) 
that maximizes the propagation length of the plasma 
waves in both mono- and multi-layer graphene struc¬ 
tures. Throughout the analysis, we consider capaci¬ 
tive effects in the device resulting from the quantum 
capacitance of graphene and the dangling bonds at the 
graphene-substrate interface. Further, apart from an¬ 
tennas and wave guides, graphene plasmons, via grated 
structures, can be profitably employed to dramatically 
increase the absorption of light at frequencies matched 
with the plasma resonant frequency. Hence, one can 
design extremely efficient terahertz photo-detectors with 
graphene heterostructures. 


Appendix A TM plasmon propagation mode in 
graphene from Maxwell’s equation 

Consider the graphene structure embedded between 
two distinct dielectric media characterized by their di¬ 
electric constants of e\ and e 2 as noted in Fig. [lO] Using 
harmonic time dependence of the electromagnetic wave 
and separating out the x —, y—, and z— components us¬ 
ing Maxwell’s equations, we arrive at the following set of 
equations. Note, only H y , E x , and E z field components 
exist as we’re considering TM propagation modes. 

\A/W— 

y 

FIG. 10: Schematic showing graphene embedded in a dielec¬ 
tric environment. The graphene layer is characterized by its 
dynamical conductivity, <j\. The incoming light is assumed 
to be p— polarized. That is, the electric held components are 
non-zero in x— and z— directions, while the magnetic held 
is non-zero in the y— direction. Surface plasmons propagate 
along x— direction. 


0E X _ . 

dz ^ 1^0 Hy 5 

(Ala) 

dHy _ . 

— lOJ€ft€ r E x , 

dz 

(Alb) 

i/3H y = —iujeoe r E z . 

(Ale) 

Further, H y must also satisfy the Helmholtz 

equation 

given as 


EE + (k 2 0 e r -p 2 )H y = 0, 

(A 2 ) 

where &o is the free-space wave-vector equal to 
the plasmon propagation constant. 

cj/c, (3 is 

For z > 0, the following components of electric and 

magnetic helds are obtained 


H y (z) = A 2 e i/3x e ~ k2Z , 

(A3a) 

rt A 

E x (z)= 2 k 2 e if3x e~ k2Z , 

a;e 0 e 2 

(A3b) 

E z (z) = ~ A2/3 e ilSx e- k2Z . 
cje 0 e 2 

(A3c) 

For z < 0, the following components of electric and 

magnetic helds are obtained 


H y (z) = A ie ilSx e klZ , 

(A4a) 












j A 

EJz) =- —kie i ^ x e klZ , 

we 0 e i 


(A4b) 


E z (z) = - —(A4c) 

cjeoei 

Now we must apply boundary conditions at z = 0. The 
tangential component of E must be matched. That is, 
E x (z = 0 + ) = E x (z = CD). Further, the discontinuity in 
the tangential component of H at z = 0 must be equal 
to J x n, where J is the electric current and n is the 
unit vector perpendicular to the surface. Doing this, we 
obtain 


A 2 — = —Ai—, 

(A5a) 

^2 ei 


. . .cr(cj,/3) . . 

A 2 A\ — i A 2 /c 2 . 

(A5b) 


cje 0 e 2 


Usin g: k 2 = f3 2 — where i = 1,2 denotes the region 
and (A5a)-(A5b), we can show that the SPP dispersion 
relation for TM waves is given as 



e 2 *c(w,/?) 



(A6) 


For the case when /3 » ^ (non-retarded regime), the 
TM plasmon dispersion relationship can be simplified to 

ei + e 2 2 icu 


{3 ~e 0 - 


r(uj,/3)' 


(A7) 


Again, separating out the real and imginary compo¬ 
nents of P and a a nd equating the real and imaginary 
components of (A7), we can show 


Pr = eo 


ei + e 2 


2ujcji 


(A8a) 



FIG. 11: Range of the validity of (A10) for TM plasmon 
propagation. 


Appendix B Electron relaxation rate in graphene 

I Scattering due to acoustic phonons 

The scattering time due to an elastic collision with 
acoustic phonon is defined as 

— = S k r(k, k')(l - cos(6>), (Bl) 

7k 

where 6 is the scattering angle between k and k' states, 
and T(k, k') is the transition rate given ad^pSl 

2 tt 

T(k, k') = — |Hk,k' \ 2 S(hvf\s! — hv fk — huj), (B2) 

where hw is the phonon energy and i7k,k' is the matrix 
element for scattering by phonons. The matrix element 
for acoustic phonon scattering is given as 


ei T e 2 / 2ujct r \ 

= (A8b) 

Using the 2D conductivity of graphene from in 
Section II, can be related to the frequency uo and the 
Fermi level Ef according to 


Pr 


2nh 2 e^e\ o 


(A9) 


where e = (ei + e 2 )/2. Re-writing the above equation 


1 l e*\E f \/3 R 

h V 27T6oe 


(A10) 


Note that the above equation is identical to ([3| in Section 
II with q standing for the real part of the propagation 
constant p R . The range of validity of the above equation 
uj/c « q < (j o/vf is demonstrated in Fig. ITT] for various 
values of e at a Fermi level of Ef = 0.4 eV. While the 
upper bound of q < u/vf is always satisfied in the plot, 
the lower bound of q » uj/c for the validity of the NR 
regime depends on the material parameters e and Ef. 


-^kjk' — COs($/2)_AqMq$k-|-q 5 k / ^ 


— iu)t 


(B3) 


where |iT q | 2 = D 2 A q 2 and |A q | 2 = h/(2pAuj cl )N(uj cl ). 
Here Da is the electron acoustic deformation potential, 
estimated to be of the order of 3£, where t is the near¬ 
est neighbor hopping integral, N( co> q ) is the occupation 
probability of phonons , an d p is the material density of 
graphene. Combining (B2) and (|B3|) with the conserva¬ 


tion of momentum (k + q = k') and inserting in (Bl), the 
final result for scattering time due to acoustic phonons is 
given as 


4 h 2 pv 2 Vf 1 

Tk = D\k^T k’ 

where v s is the velocity of acoustic phonons. In the 
above equation it is assumed that the phonon occupa¬ 
tion probability is N q ~ k R T/ Further, in the case 
of graphene, scattering due to longitudinal acoustic (LA) 
phonon modes dominates as coupling of electron-phonon 
states for other phonon modes is too weak or the energy 
scales of the (optical) phonon modes are too high for the 
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temperature range of interest. Equation (B4) can be re- kpT (see (JnJ), we obtain 


written in terms of the energy Ef according to 
4 h 3 p(v s v f ) 2 1 


Tzc — T~k F — 


D\ k B T E f ’ 


(B5) 


where kp is the Fermi wave-vector. The velocity of LA 
phonons in graphene is 7.33 x 10 3 m/s. 

II Scattering due to charged impurities 

The scattering rate due to impurities is given as 

— = A r i E k -r(k,k , )(l - cos(0 k , k ,)), (B6) 

r k 

where N{ is the number of impurities and the transition 
from state k to k' is given by the Fermi Golden rule 
according to 

r(k,k') * — I < k|V scat |k' > | 2 <5(E k - E k /). (B7) 

With the Fourier transform of the scattering potential 
V S cat (o) i the scattering rate can be written as 


h 

1~kF 


% -p—p{Ef) Ide\V scat (q)\ 2 (l-cos(0))\ (B8) 


where n\ cat is the impurities density of the scattering 
potential, q = 2kp sin(0/2), and p(Ef) is the 2D den¬ 
sity of states in graphene at the Fermi level. Using 
the Thomas-Fermi scattering potential for long-range 
Coulomb potential (see Eq. (28) nPi), the scattering 
time at Ef » kpT (large doping) is given a^M^ 


h 2 Vfk 


Emp 7~kf 


U 0 = 


fKF 


7 = 


V^fZe 2 
4e 0 e(l + 7 ) ’ 

p(£/)e 2 ? 

2eo ekp ’ 


(B9a 


(B9b) 


(B9c) 


where Ze is the net charge of the impurity atom, and nf 
denotes the concentration of charged impurities in the 
sample. Note that we have assumed the scattering to be 
elastic in nature. 

The net scattering rate is computed using the Mattheis- 
sen’s sum rule. That is, Y net = T ac + Yi mp = 


T~acEi 


of(j~ac~^~ Emp)' 


Appendix C Calculation of gate-voltage scaling law 
for maximizing plasmon propagation 
length 

Using definition of propagation length in ([5| and sim¬ 
plifying the D.C. conductivity in the case when Ef » 


Lprop — 


I e<2 \Ef\d 1 


(Cl) 


7rh 2 6oe T 

ac + I\ mp , which are computed in Appendix 


Here, T = T 

B. The relaxation rates Y ac and Y 


imp 


exhibit inverse de¬ 


pendence on the Fermi level, Ef. That is, 


Timp ~ I E f 


r ac = b\Ef\, 

where a and b are constants given as 


np Z 2 e" 


b = 


16heQe 2 ( 7rht>/eoe) 

D\k B T 


-2 


4 h 3 pv 2 v 2 ’ 


(C2a) 

(C2b) 

(C3) 

(C4) 


Substituting the definitions from ( |C2[ ) in ( |C1 ), we obtain 

(C5) 


d'prop 


/ e2 d 

E 3 / 2 

y 7iA 2 e 0 e 

a + bEj 


Noting that away from the Dirac point and neglecting 
interface trap capacitance, Ef = hvf ^itC ox /eV/, where 

C ox is the gate oxide capacitance, and V' is the differ¬ 
ence in gate voltage and the Dirac point. Using this 
relationship between Ef and U/, propagation length can 
be further simplified as 


Lprop — A 


V Q 


■/3/4 


P_V 
d v 9 J 


(C6) 


where A is a constant independent of the gate voltage, 
and (3 = b ( hvf ) 2 7reeo/e. To find the optimal value of V’ g 
that maximizes the propagation length, 


dL 


prop 


dV’ 


d 

8Vi 


v c 


■ 13/4 




= 0. 


(C7) 


Carrying out the differentiation in the above equation, 
the optimal gate voltage is given as 


v < - !£ d 
v g,opt ~ p a - 


(C8) 


Hence, one arrives at the linear scaling law between gate 
voltage and the effective oxide thickness to maximize 
the propagation length for plasmons in a gated graphene 
structure. It must be noted that parameters a and /? are 
independent of the number of layers in the multi-layer 
stack. Hence, V g opt does not depend on Ni ayer as noted 
in Section IV. 
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